Full System Engineering Design and Operation of an 
Oxygen Plant 


N93”2669^ 

/S$3 


James Colvin, Paul Schalihom, and Kumar Ramohalli 
University of Arizona/N ASA Space Engineering Research Center 



Abstract 


This paper describes one area of a project whose general aim is to produce oxygen from the 
indigenous resources on Mars. After discussing briefly the projects background and the 
experimental system design, specific experimental results of the electrolytic cell are presented. At 
the heart of the oxygen production system is a tubular solid zirconia electrolyte cell that will 
electrochemically separate oxygen from a high-temperature stream of Coleman grade carbon 
dioxide. Experimental results are discussed and certain system efficiencies are defined. The 
parameters varied include 1) the ceil operating temperature: 2) the carbon dioxide flow rate; and 
3) the voltage applied across the cell. The results confirm our theoretical expectations. 
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Introduction 


Background 

A prinrary concern for any mission to Mars must be how much energy wiii be necessary to 
complete the mission. An important consideration has to be: Shouid we continue to bring all 
propellants with us from Earth, or should we take advantage of the many known resources that are 
available to us on Mars? With this idea in mind and wishing to expand on the success of the Martian 
Viking program, Ash, Dowler, and Varsi' in the late seventies envisioned an in-situ propellant plant 
which would make use of the Martian atmosphere to produce an oxygen and methane propellant. 

The heart of this system would be an array of yttrla stabilized zirconia solid electrolyte cells. These 
cells have the ability to selectively conduct oxygen ions, thus allowing the production of pure 
oxygen. The oxygen plant has undergone many changes since it was first envisioned by Ash’ et al. 
Frisbee and Lawton have done extensive work on improving the overall system by reducing the total 
system mass and increasing the total system reliability.^'^ 

Their improved system would have the Martian atmosphere drawn in through an electrostatic dust 
Filter, which is necessary as there are numerous long-term dust storms on the Martian surface. The 
atmosphere, which consists of approximately 95% carbon dioxide, will be drawn into the system by 
a CO 2 adsorption compressor. Present pians^ require the atmosphere to be compressed from the 
ambient pressure of 6.8 mb to a pressure of 1 bar for delivery to the cathode of the electrolyte. 
Before entering the zirconia array, the flow will pass through a heat exchanger which will raise the 
temperature from the COj compressor's exit temperature of 600 K to approximately 1000 K. The 
source of the heat for the heat exchanger is the exhaust flow from the array. Once the flow has 
entered the array, it will be further heated to a temperature of 1273 K. 

This temperature is sufficient to partially dissociate the carbon dioxide into carbon monoxide and 
diatomic oxygen. The O 2 will be increased to 4 bar. The cell's exhaust will consist of mostly CO 2 
with some CO resulting from the removal of O 2 . The O 2 produced at the anode will pass through 
a radiator where it will be cooled from 1270 to 250 K with a pressure of 3.8 bar. The O 2 will next 
pass through an O 2 adsorption compressor where its pressure will be Increased to 28 bar and its 
temperature to 400 K. After passing through a radiator, the flow will be cooled to 230 K. The O 2 will 
be finally cooled to 100 K by a molecular adsorption cryo-cooler refrigerator and stored for its final 
use. The use could initially be the oxidizer for the propellant necessary to return a Martian sample 
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to earth, and then eventually, life support for a manned Martian mission. 


Yttrla-Stabillzed Zlrconla Solid Electrolyte 

In Lawton's work,^ he lists development risk Actors for components. In his option III, the oxygen 
cell is the only component still listed as risk ^ctor 4 , meaning Ihere are still serious problems that 
must be addressed as well as some intensive development required.’ This is an area of current 
research. Further details were worked out in Reference 5 . Since the Martian atmosphere is 
predominantly COj, the remainder of this report will refer to the atmosphere as COj. The COj, when 
it enters the cell array, will be heated to 1273 K At this temperature, the CO2 will begin to partially 
dissociate into CO and Oj. The zirconia electrolyte is sandwiched between two porous platinum 
electrodes. The dissociation will occur at the cathode with the O2 entering the electrode and moving 
towards the electrode-electrolyte interface. The driving force for this movement is the partial 
pressure gradient developed by the electrolyte removing oxygen from the interface area. Once the 
O2 reaches the interface, it is further reduced to monatomic oxygen. The oxygen received two 
electrons from the negative electrode and becomes an oxygen ion and begins to migrate through 
the zirconia electrolyte towards the anode. Upon reaching the anode, the ion will release its two 
electrons to the positive electrode and then recombine with another oxygen atom to reform the O2 
molecule. Qualitatively, this describes the oxygen separation process. Figure 1 * shows this process 
schematically. 

Richter performed the initial intensive testing of this el^rochemical process with the aim of 
quantifying this physical procedure.^ His work was performed using a tubular zirconia ceil. He 
developed the basic thermodynamic and electrochemical models for the reduction of CO2 and the 
subsequent production of O2. A few years later. Suitor continued the investigation.^ In his 
experiment, he used the disk geometry for his cells. Additionally, he investigated the use of different 
electrode materials. Although his supply gas was air, many of his results can also be applied when 
CO2 is the supply gas. 

In this study, zirconia cells with the tubular geometry were used. The aim was to develop various 
cell efficiencies by varying several system control parameters. These parameters were: the potential 
applied across the electrolyte, the electrolyte operating temperature, and the incoming CO2 flow 
rate. We would like to know how the CO2 production rate varies as a function of these parameters. 
The various efficiencies include the cell's Nemst efficiency (a measure of the theoretical energy 
required against the actual energy put into the system). The definitions of these efficiencies will be 
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C02 -> CO + 02 ( Pi ) 



02 ( P2 ) 

Fig. 1 Oxygen migration through the ZrOj cell. 



Thermocouple Pressure Transducer 

© ® ® Flow/VoU/Current Meter 

Fig. 2 Single cell test bed (schematic). 
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developed later in this report. The singie zirconia ceil test bed currently being studied is shown 
schematically in Figure 2. The test bed consists of the following: 1) one tubular zirconia celi; 2) one 
voltage DC power supply; 3) two digital multimeters; 4) two ceramic dam shell heaters; 5) one 
Watlow heater controller; 6} Kaowool type ceramic Insulation; 7) one Oj pressure transducer (0-500 
psia); 8) one CO 2 pressure transducer (0-500 psia); 9) one O 2 flow meter (0-50 seem); 10) one COj 
flow meter (0-5000 seem); 11) three ktype thermocouples; (12) Coleman grade CO 2 for supply gas; 
13) one PC/386 processor for data acquisition; and 14) Varian model 3700 gas chromatograph. 

To simulate the Martian atmosphere, which contains approximately 95% CO 2 and only 0.13% O 2 , 
Coleman grade CO 2 was used to supply the test bed. This grade contains less than 20 ppm O 2 . The 
flow was maintained at slightly above local atmospheric pressure (13.7 psia) at a temperature of 
75°F. The flow rates were varied between 38 and 1475 seem. Row would enter the zirconia tube 
through a 1/8 in. alumina tube and then pass to the far end of the zirconia chamber (see Figure 3). 
The clam shell heaters are centered about the middle 7 in. of the zirconia device. This means 
hearing of the flow will begin in the alumina tube. Exiting the alumina tube, the flow reverses 
direction while continuing to be heated and fiov^ across the cathode of the electrolyte. The CO 2 , 
now heated, begins dissociation and is drawn to the cathode. The free stream is now a mixture of 
CO 2 , CO, and Oj. Oj is dissociated to monatomic oxygen at the electrode-electrolyte interface and 
eiectrochemically conducted through the zirconia to the anode while the COj and CO exhaust pass 
out through the exit of the tube. The supply CO 2 flow rate is controlled by a metering valve in the 
exhaust flow. The O 2 produced (at approximately 13.7 psia) flow« through the mass flow meter and 
then can be directed to a water bubbling device, a sample cylinder, or directed to the gas 
chromatograph for analyzing. The zirconia electrolyte and its platinum electrode have an upper 
temperature limit of 1 150°C where a phase change in the zirconia will take place. A critical voltage 
limit of 2.23 VDC was assumed in accordance with the work of Frisbee.^ Keeping below this 
F>otential will prevent the oxygen from being driven from the zirconia lattice structure causing 
permanent damage to the cell. With these material limits in mind, a self-imposed limit of 1 100°C, and 
2.0 VDC was used during all testing. 

Single-Cell Test Results 

Fundamental Results 

A series of tests were conducted to attempt to characterize the effects of temperature, cell 
potential, and carbon dioxide flow rate on the production rate of oxygen. During testing, the 
temperature was varied from 800 to 1100°C in increments of 25 q. Cell potential varied from 
approximately 0.6 to 2.0 VDC In increments of 0.1 volts. The CO 2 flow rate varied from 38 to 1475 
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seem. The fdlowing figures show the results this testing. Figure 4 shows the dependenee of the 
produetion rate of Oj on the applied cell potential for a variety of COj flow rates at a temperature 
of 1000°C, while Figure 5 displays the production at four different temperatures at a CO 2 flow rate 
of 138 seem. Note, on both Figs. 4 and 5, the second-order dependence of oxygen production on 
the applied cell potential. 

Figure 6 demonstrates the dependence the oxygen production has on the temperature of the cell 
with an applied cell potential of 2.0 VDC. Notice the weak dependence the oxygen production has 
on the carbon dioxide flow rate. Figure 7 presents the oxygen production vs cell temperature at a 
carbon dioxide flow rate of 138 seem. Figure 7 also dearly depicts the dependence on voltage of 
the oxygen production rate. 

Figure 8 indicates the O 2 production rate dependence on the CO 2 supply flow rate for an applied 
cell potential of 2.0 VDC. This is the most graphic Illustration of the lack of dependence on the 
carbon dioxide flow rate for the oxygen production rate, especially for CO 2 flow rates greater than 
about 200 seem. Figure 8, like the previous figures, shows the dependence on the temperature for 
the O 2 production rate. 

Interpretation of the Results 

The results of the extensive system analysis can be discussed in five areas; 1) flow rate effects, 2) 
mass flow ratio, 3) oxygen conversion efficiency, 4) Nemst efficiency, and 5) the system efficiency. 
Each of these areas will be addressed in order. 

The flow rate effects were the easiest to interpret. The basic schematic is shown in Figure 2. It was 
initially found that the oxygen yield rate (production rate) increased as the CO 2 flow rate increased, 
but only up to a certain point, after which increased flow rate actually resulted in the decrease of 
the yi^d of oxygen. This was observed while other parameters such as cell voltage and cell 
temperature were held constant. It was suspected that the carbon dioxide may not have had 
sufficient residence time within the heated cell to achieve the required temperature at the high flow 
rates. To verify this hypothesis, a simple thermal diffusivity analysis was performed, after confirming 
that the flow within the tube Is indeed laminar via simple Reynolds number calculation. The 
residence time is given by i/u, where / is the tube length and u is the mass averaged velocity. The 
characteristic time for heat transfer within the tube is given by (cf/2)^/a, where d is the cell inside 
diameter and oc is the thermal diffusivity {k/pc^ of the carbon dioxide at high temperatures. As can 
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be seen in Figure 9, the peak occurs almost precisely where the two times become equal. It is thus 
dear that up to the point where the residence time Is sufficient to bring the carbon dioxide to the 
requisite temperature, the production rate of the oxygen increases with the flow rate, but beyond 
this point the production rate fells off. It would be instructive to measure the actual gas temperature, 
and ncA merely the cell wall temperature as Is done in the experiments reported here. 

The mass flow ratio (the mass production rate of Oj divided by the COj supply mass flow rate) is 
of much interest. This is due to sizing constraints within scaled-up systems such as the primary 
carbon dioxide compressor, the heat exchanger, and the radiators. For these reasons. Figs. 10 and 
1 1 plot the mass flow ratio vs COj flow rate and temperature, respectively. 

In previous studies of the propellant production plant, the emphasis was placed on sizing of the total 
production system. For these studies, oxygen conversion efficiencies of 25-30% were used when 
considering the estimated surface area needed to produce a required amount of oxygen per day. 
The use of these efficiencies was necessary to minimize oxygen plant mass. However, when 
considering the operation of the zirconia cell alone. Figure 12 shows the complete range of oxygen 
conversion efficiencies. 


As the concern of this report is the complete characterizaticHi of the zirconia electrolytic cell, all 
conversion efficiencies must be investigated. Richter, in his early work® on the reduction of CO 2 
began with the basic Nemst relation; 


— ln[ — 
^ zF X, 


02 


Pi 


Here, is the Nemst voltage, T is the temperature (K), z is charge/mole (z = 4), /? is the universal 
gas constant, and F is the Faraday constant. Additionally, Pj is the pressure at the anode while 
P, is the partial pressure of the oxygen In the supply COj flow stream at the cathode. Richter then 
continued with his model and developed an instantaneous Nemst potentiefl 


E„,-jgln[ ■ ./s -; 1 


zF ^ K{l-n) ^2 


n 


IV-41 





IV-42 



EXPLANATION OF FLOW-R ATE EFFECT 



[a is th* tharmal diffusivlty of CO, at 1000* C] 


Fig. 9 Thermal diflhsivity vs residence time. 



Fig. 10 Oj/CO, mass flow ratio vs CO, flow rate, cell potential 
2.0 VDC. 
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In this equation, K is the normal equilibrium constant and n is the mole fraction of COj reduced to 
CO. The value n can be calculated from this equation: 


n-C{ 



The value C is a constant from Richter's work (C = 6.969), / is the measured current produced 
across the ceil, and Vco 2 is the volumetric flow rate of the supply CO 2 stream. The current (A) was 
measured experimentally, as well as the volumetric flow rate of the CO 2 (seem). The cell 
temperature (K) was measured by a k-type thermocouple positioned at the outside wall of the cell, 
and the pressure of the O 2 collection side of the cell was assumed to be one atmosphere. Using 
the proper values for the equilibrium constant K.® there was enough information available to 
calculate the instantaneous Nernst voltage. Richter cautions the use of this equation by stating that 
it is veilld only if the critical voltage over the cell is not reached. The critical voltage can be defined 
as the potential which will be just sufficient to begin removing oxygen from the lattice of the zirconia. 
This process would be apparent from an "elbow’ being observed in the current vs potential plot in 
the higher potential range (> 2 VDC). Figure 4 shows this potential has not been reached due to 
the lack of an elbow. According to Richter then, the potential can be expressed as a function of the 
operating temperature and the actual amount of O 2 removed and not the partial pressure of the O 2 
at the cathode. The instantaneous Nernst potential was used in defining the Nernst efficiency; 


This equation shows the Nernst efficiency is the ratio of the theoretical potential divided by the 
actual potential {E^J. What this indicates is the instantaneous amount of cell overpotential. Figure 
13 shows the results of this analysis. 

It is apparent from this figure that the amount of cell overpotential rises as the applied cell potential 
rises. One source of this overpotential could be the pressure drop through the porous electrode. 
It is important to remember that from Figure 4, the rise in O 2 production is almost directly 
proportional to the potential applied to the cell. Ash® et al. point out that experiments have shown 
that as the O 2 production rate rises, there is an Increase in the pressure drop across the negative 
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electrode which can approach the partial pressure of the oxygen In the CX )2 stream resulting in an 
ever increasing pumping power requirement. Another possible source of the overpotential is slow 
diffusion, adsorption, or dissociation processes near the electrode-electrolyte interface of the cell 
resulting in concentration overpotentlal.^'’ This report also states the possibility of a transition 
overpotential due to slow electrochemical reaction. The sources of overpotential cannot be 
fficperimentally determined with this system as currently constructed, but any future characterization 
of the zirconia electrolyte should definitely seek to identify the source of these overpotentials. 

With the Nernst efficiency established, it was next desired to arrive at a system efficiency. The 
contrd volume is taken around the entire cell system. The useful work inside the control volume Is 
divided into two parts. The first being the rate of thermodynamic energy being used to dissociate 
the CO 2 . The second being the work required to electrochemically conduct the oxygen ions through 
the zirconia. The power passing into the control volume is also divided into two components. The 
first being the power delivered to the ceramic heaters, arxJ the secorxJ being the power delivered 
to the electrolyte. The energy required to thermally dissociate the CO 2 is given by the equation: 

p^l ■‘■■^ 7’“ -^29 8 ^ ^ j 


Here, Qjh Is the thermodynamic power input with the calculation showing the normal enthalpy of 
reaction equation with the enthalpy of formation and sensible enthalpies. The subscripts / and / 
represent the products and reactants respectively. The coefficient N represents the molar flow rate 
of each respective constituent. This value and the other power values give the following definition 
for the system efficiency: 

_ TH 

^ sy 8 rp Y 4 - S’ T 

^ cell ^ heaf^ heat 


In this equation, the subscript heat refers to observed values of the ceramic heater, and the 
subscript cell refers to observed values with the electrolyte, E being the potential (V) and / being the 
current (A). Figure 14 displays the results of this analysis. 

Summary 

This paper has presented results from the first phase of a well-planned multiphase research program 
aimed at significantly reducing costs of future space missions; terrestrial applications are always 
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kspt in nnind. This proof of concept study has used full system hardware, realistic solid electrolytic 
cell material, realistic operating temperatures, full cell voltages, and full scale flow rates. These 
features distinguish our experiments from the usually understood test-tube demonstrations, where 
the conditions are substantially different from the real life counterpart. The only feature that is not 
duplicated here is the scale of oxygen production. 

Important efficiencies are defined and measured. The basic electrolytic efficiency refers to the 
efficiency of using electrical input in producing oxygen, while the overall system efficiency refers to 
the enthalpy difference achieved between the product stream and the reactant stream divided by 
the overall energy input. The carbon dioxide flow rate, cell voltage, and cell temperature are all 
varied parametrically. A simple heat transfer theory explains the flow rate effects. The voltage effects 
are in good agreement with the manufacturer's specifications. 

These thorough characterizations of the component performance parameters naturally lead us to 
the next step of scale-up and full system demonstration with active controls. Creative solutions to 
these engineering designs indicate that future space missions could realize substantial cost savings 
through the use of local (in-space) resources. 

Acknowledgments 

Research was supported by NASA Code R, through the University of Arizona Grant NAGW-1332. 
The authors are grateful to Murray Hirschbein and Gordon Johnston tor funding and helpful 
comments. 


References 

^Ash, R., Dowler, W.. and Varsi, G., Feasibility of Rocket Propellant Production on Mars," Acta 
Astronautica, Vol. 5, 1978, pp. 705-724. 

^Lawton, A.L., and Frisbee, R.H., "A New Look at Oxygen Production on Mars ISPP," Jet Propulsion 
Lab. Paper D-2661, Sept. 1986. 

^Lawton, A.L, "Risk Factors in the Development of Zirconia Cell Technology for the Production of 
Oxygen from the Martian Atmosphere," Jet Propulsion Lab. Paper D-3546, Aug. 1986. 

“Frisbee, R.H., "Mass and Power Estimates for Martian In-Situ Propellant Production Systems," Jet 
Propulsion Lab. Paper D-3648, Oct. 1986. 

^Ramohalli, K., Lawton, E., and Ash, R., "Recent Concepts in Mission to Mars: Extraterrestrial 
Processes," Journal of Propulsion and Power, Vol. 5, No. 2, 1989, pp. 181-187. 


IV-48 



®Richter, R., “Basic Investigation into the Production of Oxygen in a Solid Bectrolyte,’ AIAA Paper 
81-1175, June 1981. 

^Suitor, J.W., Berdahl, C.M., Ferredl, J.F., Marner, W. 

Schroder, J.E., and Shichta, P.J., "Development of an Alternate Oxygen Production Source Using 
A Zirconla Solid Electrolyte Membrane," Jet Propulsion Lab. Paper D-4320, May 1987. 

®Wark, K., Thermodynmnics, 3rd ed., McGraw-Hill, New York, 1977, p. 840. 

^Ash, R., Richter, R., Dowler, J.A., Hanson, and Uphc^, C.W., "Autonomous Oxygen for a Mars 
Return Vehicle," International Astronautical Federation Paper 82-210, 1982. 

’°Etsell, T.H., and Flengas, S.N., "Overpotential Behavior of Stabilized Zirconia Fuel Cells, Journal 
of die Electrochemical Society, Vol. 118, No. 12, 1971, pp. 1890-1900. 


IV-49 



